A=[0 0.5 0 0.5;
   0.5 0 0.5 0;
   0 0.5 0 0.5;
   0.5 0 0.5 0];
disp('二次型矩阵为：');
disp(A);
[P,T]=schur(A);
disp('正交变换矩阵为：');
disp(P);
disp('二次型的标准型为：');
disp(T);